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The magnetotransport of two dimensional holes in a GaAs/AlGaAs heterostructure is studied 
experimentally and theoretically. Spin-orbit splitting of the heavy hole band is manifested at high 
carrier densities in two Shubnikov-de Haas frequencies, classical positive magnetoresistance, and 
weak antilocalization. The latter effect combined with inelastic scattering between the spin-orbit 
split bands lead to metallic characteristics, namely resistance increase with temperature. At lower 
densities, when splitting is smaller than the inverse elastic scattering time, the two bands effectively 
merge to yield the expected insulating characteristics and negative magnetoresistance due to weak 
localization and interaction corrections to the conductivity. The "metal to insulator" transition 
at intermediate densities is found to be a smooth crossover between the two regimes rather than 
a quantum phase transition. Two band calculations of conventional interference and interaction 
effects account well for the data in the whole parameter range. 
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I. INTRODUCTION 

The study of magnetotransport in two dimensional 
(2D) electronic systems regained considerable interest 
since the observation of Kravchenko et alB that the re- 
sistance of high mobility, high density silicon MOSFETs 
decreases and saturates to a residual value as the tem- 
perature is reduced. This metallic characteristic is in 
sharp contrast with- the prevailing dogma that 2D sys- 
tems are insulating3 although interaction may possibly 
lead to delocalizationfl The conflict between Ref. |^ and 
the expected insulating characteristics motivated exten- 
sive experimental and theoretical efforts. Soon, simi- 
lar characteristics have been flbserved in other silicon 
samples,El SiGe quantum wellsja AlAs based 2D electrop, 
gas (2DEG),Q 2DEG in GaAs/AlGaAs heterostructures, 
2DEG in GaAs with self assembled InAs quantum dots, 
and various realizations of 2Ijl||4jolft— gas (2DHG) in 
GaAs / AlGaAs heterostructures £lll3'Bt3 

The samples that show metallic characteristics share 
some other features: 

(a) Insulating behavior at low carrier densities, namely, 
resistance increase with decreasing temperature. 

(b) Weak dependance of the resistance upon temperature 
at some intermediate density (coined "metal-insulator 
transition" (MIT)). 

(c) Negative magnetoresistance at the density corre- 
sponding to the MIT. 

(d) Large, positive magnetoresistance for magnetic fields 
parallel to the layer. 

(e) A crossover from metallic to insulating characteristics 
for large enough parallel magnetic fields. 

The zero field crossover from metallic to insulating 
behavior was identifie d . by . [-»3mc authors as a second 
order phase-jtransitionJldllil'La Scaling theory has been 
constructedtj and even reentrant transition to an insu- 
lating phase, at still higher densities, has been argued to 



occur for holes in GaAs.0 

Notwithstanding the remarkable similarities between 
the magnetoresistance and temperature/density depen- 
dences of the different material systems and samples 
listed above, some gross differences should be appreci- 
ated: 

(i) While the weak field magnetoresistance in silicon sam- 
ples is negative in all regimes, it is positive for holes in 
GaAs in most of the metallic regime (high carrier con- 
centration) . 

(ii) While only single band is observed in silicon sam- 
ples, the 2DIIG in GaAs samples, in most of the metallic 
regime, is characterized by two distinct bands. These 
bands are manifested in two Shubnikov-de Haas (SdH) 
frequencies, as well as in classical, positive magnetoresis- 
tance. 

(iii) The relative resistance change with temperature in 
different materials vary between about 1% for electrons 
in GaAs to an order of magnitude in silicon MOSFETs. 

(iv) The "critical resistance" at the MIT varies substan- 
tially from sample to sample and between different ma- 
terial systems. In contrast to some early claims iL can 
deviate substantially from one resistance quantumti., 

(v) The large effect of back gatirw, in GaAs 2DHG0'EI 
is abserd_for electrons in silicones or is at least very 
differentB. 

The extent to which the magnetotransport in different 
2D systems is universal is hence unclear at this point. 
The multivalley band structure of silicon, the short range 
potential fluctuations in silicon as opposed to the long 
range ones in GaAs, the strong spin-orbit coupling in 
III-V materials, not to mention obvious differences in ef- 
fective masses and Zeeman factors, might all turn out to 
be important. 

A wide spectrum of mechanisms has been proposed 
to account for these-.features, including a new type 
of superconductivityJ13 a novel metallic phase induced 



2 



by disorder enhanced interactionJli (this direction is 
based on the earlier jwprks in Ref. [sj), Wigper glass 
and non-Fermi— ]iquid,t3 a new liquid phaseJ23 irnpn-, 
rity scatteriaaS temperature dependpit screening,Efl'E£l 
spin efj[e£±s,EirE3 interband scattering, band structure 
effects,tlrE3 and classical percolationE] None of these ex- 
planations account for all features in all materials. It 
would thus be fair to state that the metallic behavior, 
as well as the other features listed above, remain unex- 
plained as universal phenomena. 

The present manuscript presents an extensive experi- 
mental and theoretical study of magnetotransport in two 
dimensional hole gas in GaAs/GaAlAs heterostructures. 
The experimental data pertain to measurements of high 
mobility samples as well as low mobility ones. The data 
display all the characteristic features listed above. At 
high densities all samples display metallic characteristics 
(Figs. ^, ||c-e). At low densities the samples are insulat- 
ing (Fig. 6a) in accordance with point (a) above (we have 
chosen to display the insulating data for the low mobil- 
ity sample to avoid possible inhomogeneities that might 
occur at the very low densities where the high mobil- 
ity samples turn insulating). At an intermediate density 
(Fig. Hb), the samples display weak dependence of the 
resistance upon temperature ("MIT", point (b) above). 
All samples are characterized by a large, positive magne- 
toresistance for a parallel magnetic field (points (d), (e) 
above and Fig. |^) and suppression of the metallic char- 
acteristics by such a field. 

In a recent letteiEl we reported magnetotransport 
measurements done on a high mobility 2DIIG sam- 
ple. Concentrating on the high density regime we were 
able to show that the metallic characteristics result 
from inelastic scattering between the two heavy hole 
bands split by the spin-orbit interaction. The posi- 
tive magnetoresistance observed in that sample was fully 
accounted^[foiL-.by the well known classical two band 
formulaecSEBo and the resistance increase with tem- 
perature was shown to result from the enhancement 
of Coulomb scattering between the two spin-orbit split 
bands. That work thus explained the metallic character- 
istics with well known semiconductor physics without in- 
voking a novel metallic phase or any other "new" physics. 
The spin-orbit ac]i±|haiids,were studied extensively. -bpth 
experimentallyjBll|£ilE3Ej and theoretically.oHE30 In 
fact, Murzin et alSB have previously shown that inter- 
band Coulomb scattering is manifested in a resistance 
increase with temperature. 

The present manuscript extends our previous work 
in several ways. First we refine the high density data 
analysis and show that the negative magnetoresistance 
at larger magnetic fields results frorn miantum correc- 
tions due to Coulomb interaction.ESEJo Next we re- 
port new measurements covering the insulating, MIT, 
and metallic (high density) regimes. The data con- 
form W|it h || pr j OV | io iL is . piag petotransport measurements on 
2DHG.BByyka.BBEl Careful analysis then shows 
that well-known physics accounts for all data. The 



emerging picture is briefly as follows 

It is well 
experimental! 



theoreticallyilEiil'i and 
that the bulk heavy hole 
band in assymmetric GaAs/AlGaAs heterostructures is 
split by spin-orbit interaction. The resulting sub-bands 
are manifested by two SdH frequencies, as well as by a 
classical, Lorentzian shaped, positive magnetoresistance 



(PMR). The band splitting energy, eg, vanishes at k = 
and increases with fco At low densities, when eg at the 
Fermi energy is small and the elastic scattering time, r, 
is short, eg <C fi/r, the two bands are strongly mixed 
and spin-orbit splitting is unimportant. The resulting 
magnetoresistance is negative and dominated by weak lo- 
calization (WL) and hole-hole Coulomb interaction. As 
the temperature is increased, these quantum corrections 
are suppressed and the sample becomes more conductive. 
At high densities, the sub-band splitting is large and t is 
long, thus eg ^ h/r. The two sub-bands are hence dis- 
tinguishable. The weak field magnetoresistance is pos- 
itive and dominated by weak anti-localization (WAL) 
induced by inter sub-band scattering (effectively, spin- 
orbit scattering), hole- hole interaction (Altshuler-Aronov 
correctionl2j) , and classical two band magnetoresistance. 
A slightly higher magnetic field suppresses the quantum 
corrections in the Cooperon channel while the classical 
magnetoresistance persists. At still higher fields, the clas- 
sical magnetoresistance is also suppressed and the mag- 
netoresistance turns negative due to hole-hole interaction 
in the diffuson channel. 

The metallic characteristic, dp/dT > 0, in the high 
density regime is traced to the enhancement of inelastic 
inter sub-band scattering as well as the suppression of 
quantum corrections (especially the WAL) with increased 
temperature. The weight of quantum corrections relative 
to the classical effect depends on sample resistivity. The 
quantum corrections are negligible in high mobility sam- 
ples and very important in the low mobility ones close 
to the MIT. The weak temperature dependence of the 
resistivity at the MIT merely reflects effective cancella- 
tion of the resistance increase with temperature due to 
inelastic interband scattering and suppression of WAL, 
against resistance decrease with temperature due to the 
suppression of hole-hole interaction. 

Our analysis thus explains how the sample crosses 
over from an insulator at low densities to a "metal" at 
higher densities. Moreover, we account for the change 
from negative to positive magnetoresistance as the den- 
sity is increased. All these features emerge from the well 
known band structure of holes combined with conven- 
tional quantum corrections to the magnetoresistance. No 
new metallic phase is invoked. 

The manuscript comprises the following parts. Section 
II pcoyides a refined analysis of our previously published 
datall3 on high mobility samples. The refined analysis in- 
cludes hole-hole interaction in the diffuson channel. The 
extracted inter and intraband scattering rates remain un- 
changed. The analysis provides new information about 
the Coulomb interaction parameter at different densities. 
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FIG. 1: Hole density in the light (pi) and heavy (ph) spin- 
orbit split bands as a function of total density. Inset - One of 
the Shubnikov de Haas traces (high mobility sample) used to 
determine the hole densities. 
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FIG. 2: The two-band classical contribution to the magne- 
toresistance for different temperatures (high mobility sample). 
Note the perfect Lorentzian shape MR at low fields. Inset- 
Full magnetoresistance curves for the same temperatures. 



The theory of quantum corrections to t he magnetotrans- 
port is briefly reviewed in Section II E . Previous calcu- 
lations are extended to account for two sub-bands with 
different mobilities. Details of these calculations are pro- 
vided in the Appendix. Extensive data measured on a 
lower mobility sample are provided in section III. The 
data are analyzed taking into account WL, WAL, hole- 
hole interaction in the cooperon and diffuson channels, 
and interband scattering. A comprehensive picture of 



magnetotransport in a 2DHG is compiled in section IV. 
Section ^ details some theoretical aspects of magneto- 
transport in a two band system with spin-orbit splitting. 
Plasmon mediated Coulomb scattering is analyzed in this 
section, and shown to yield an Arrhenius temperature de- 
pendence. Finally, this section also discusses the quan- 
tum corrections to the magneto-conductivity in this sys- 
tem. 



II. HIGH MOBILITY 2DHG SAMPLE AND 
THEORETICAL BACKGROUND 



spectively, not to be confused with the bulk light and 
heavy bands). Fig. ^ depicts pi and ph as a function of 
the total density, Ptotah for the high (circles) and low 
(squares) mobility samples. The band's population in 
the two samples are similar. Below a total density of 
about 2 X 10^^ cm~^ it is hard to resolve two bands in 
the SdH data. For higher densities, the bands split and 
for ptotai > 4 X 10^^ cm~^ practically all additional car- 
riers populate the heavier, less mobile, band. 

The weak field MR of the high mobility sample at 
Ptotai = 4.25 X 10-" cm-2, pi = 1.52 x 10"" cm'^, 
Ph = 2.73 X 10^^^ cm^^, and different temperatures is 
presented in the inset to Fig. ||. The Lorentzian shaped 
PMR expected from two band transport is obtained by 
subtracting the predicted quantum correction (see below) 
from the full MR data. The result is depicted in Fig. |[ 



A. Experiment 

The high mobility sample had a mobility of ~ 
300,000 cmV(Vs) at p = 4 x 10" cm^^ and T = 
200 mK. The low mobility sample had a ~ 
20, 000 cm^/(Vs) mobility at a similar carrier density 
and T = 400 mK. The 2DHG was confined in both sam- 
ples to a GaAs/Alo.8Gao.2As interface in the (100) plane. 
The samples had a 2DEG front gate, 40 and 20 nm (high 
and low mobility samples, respectively) above tho. 2DHG 
and a silicon doped back gate, 300 nm below iiS3. 

The inset to Fig. |^ depicts a characteristic SdH curve. 
Two distinct frequencies corresponding to two bands are 
clearly observed and used to determine the carrier den- 
sities in the lighter and heavier bands {pi and ph, re- 



The analysis of the high mobility data presented 

ilia 

was based on standard, two-band transport fo™jiulae,E2 
generalized to include interband scattering. c2lE2l The 
Lorentzian PMR expected for two band transport was 
extracted by subtracting the weak parabolic nega- 
tive magnetoresistance (NMR), attributed to Coulomb 
interactionsE3, from the full MR curve. Here we refine 
the analysis and incorporate quantum corrections ab ini- 
tio. These corrections include hole-hole interaction, WL, 
and WAL. They all become increasingly important as 
the sample conductivity is reduced. Their inclusion ab 
initio facilitates a unified analysis of high and low mobil- 
ity samples. Moreover, the quantum corrections provide 
valuable information on the fundamental scattering pro- 
cesses in 2DHG. 
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B. Theoretical background 

Our analysis generalizes Choi et a/.0 approach to the 
case of two band transport. The starting point is a 4 x 4 
conductivity matrix for the two bands in the presence 
of a perpendicular magnetic field. All matrix elements, 
except the Hall ones, are modified by quantum effects. 
Setting the Hall currents to zero and inverting the matrix 
one obtains the longitudinal resistance with WL, WAL, 
and interaction corrections: 

pUH) = PLiH) + fiSi,Sh,Q,H)-{5a,,-^+Sa„u), (1) 
where 

L 



1 + {H/Wf 



Pl{H oo) = 



W = 



L 



R^Si + RfSh — 2QRiRh 

{Ri + Rh? ' 
Si + Sh + 2Q 
Ri + Rh 

[Ri{Sh + Q)-Rh{Si + Qt 
{Si + Sh + 2Q){Ri+Rh)^ ' 



(2) 
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Here, l^h correspond to the lighter and heavier heavy 
hole sub-bands, respectively. Ri — 1/epi, with i = l,h, 
is the Hall coefficient of the i-th band. The diagonal 
resistances Si, Sh, and off-diagonal resistance Q can 
be expressed in terms of elastic and inelastic contri- 
butions, SiiT) = SiiO) + a"i[Q(r) - Q(0)]; ShiT) = 
Sh{0) + a[Q{T) — (5(0)], where a is a function of ve- 
locities and densities. The resistances 5/(0), Sh{0) per- 
tain to inter and intraband impurity scattering, Q{0) 
and Q{T) — Q{0) pertain to elastic and inelastic scat- 
tering, respectively. The latter processes may include 
carrier transfer between bands as well as drag-like pro- 
cesses where a particle from one band scatters off a par- 
ticle in the other band and both carriers maintain their 
bands. The resistance p^, reflects classical, two band 
PMR and depends in a Lorentzian way on the magnetic 
field. The conductances, Sa^L and Saint, stand for WL 
(or WAL) and interaction corrections to the conductiv- 
ity. The function f{Si,Sh,Q,H) given in the appendix, 
is roughly quadratic in H. It is negative for H < W and 
positive for H > W (see Fig. |5|). The functions Sa^L, 
Saint, and f{Si,Sh,Q,H) depend on both temperature 
and the magnetic field. Equation (Q) shows that the over- 
all quantum contribution to the resistance comprises the 
corrections to the conductivity of both bands. A qual- 
itative understanding of / can be gained from the sin- 
gle band case where it reduces to — [l — (^c^)^] , with 
ujc = eH/mc being the cyclotron frequency [cf. Eq. (Ju)]. 



The quadratic increase of / with H yields a large inter- 
action correction to the magneto-resistance despite the 
smallness of Saint- 

We turn now to review the theory for the quantum 
corrections, WL, WAL, and interaction, taken into ac- 
count in the data analysis. Further details are pro- 
vided in Scc^j.^ J¥e start with WL and WAL. Vari- 
ous author£30tlj'LZl calculated the hole band structure 
of asymmetric GaAs/AlGaAs quantum wells and het- 
erostructures and found that each band is characterized 
by a wave-vector dependent spinor state. When a hole 
is scattered from k to k' the spin-orbit couolijpg, leads 
to Dyakonov-Perel spin precession around k' L3c3't3'OC3. 
As a result, the spin relaxes at a rate Ti/tso ~ e^r/h, 
where r is the momentum relaxation time, and eg is the 
energy difference between the two bands for an average 
kp. Since WAL is important for g ~ 1, where the momen- 
tum relaxation rate is similar to the single particle scat- 
tering time, we neglect the difference between the two. 
Scattering between two spin-orbit split bands is hence 
equivalent to spin-orbit scattering by, e.g., a large atom. 
For Tso < {t^ is the dephasing time), WAL takes place 
with either positive or negative MR for = ^JficjeH 
larger or smaller than Isofj \/Dtso {D is the diffusion co- 
efficient). Averkiev et alxA calculated the WL correction 
for weak and strong spin-orbit interaction in symmetric 
p-type quantum wells. They found that the key param- 
eter for the anomalous MR is kpa/ir, where a is the well 
width. This parameter is a measure of bulk heavy/light 
hole wave function mixing. For kpa/n ~ f they predict 
PMRr-at weak fields and NMR at higher fields. Pedersen 
et alsB measured such effe ct in a symmetric GaAs p-type 
quantum well. In section |VB| we generalize Averkiev et 
al. procedure to the case of asymmetric GaAs p-type 
quantum wells with their spin-orbit split bands, taking 
into account interband and intraband particle-particle 
propagator contributions to the Cooperon equation. The 
final result is 



Sa,,AH) = ^wL{H)-aoiH) 



(3) 
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where "if is the Digamma function, and the characteristic 
magnetic fields Htr, Hg^ and H^p are given by 



Htr — 



he 
AeDr 



Hso — 



Tic 



-ieDTg. 



eir/h, it is 



Using Dyakonov-Perel equation, h/Tso 
then possible to estimate the spin-orbit band splitting, 
tg, from transport measurements. The existence of two 
bands is manifested in three phenomena; a Lorentzian 
PMR, two SdH frequencies, and finally, WAL correction 
to the conductivity. While the SdH frequencies merge 
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into a single frequency when the gap becomes smaller 
than the scattering rate (e^ < ?i/t), the WAL remains 
visible down to considerably smaller gaps provided the 
temperature is low enough. The WAL is hence a power- 
ful tool in the analysis of the 2DHG band structure and 
scattering mechanisms. The WAL is visible as long as 
Tso < T^- Substituting h/r^ —T/g for the dephasing rate 

we find WAL should be visible as long asT < (~|~)^' 
although tgTjfi may be smaller than unity. Moreover, 
since the relevant magnetic field for WAL is g times 
smaller than that relevant for the classical effect or the 
SdH oscillations, the WAL probes the_Lwo bands at very 
small magnetic fields. Winkler et a/.EJ argue that SdH 
data at finite fields underestimates the zero field band 
splitting. Our analysis proposes WAL as a better tool 
for studying degeneracy lifting of the heavy holes, par- 
ticularly at weak magnetic fields. The WL correction 
is significant as long as Ih > I or \/2u>cTg < f, where 
ujc = eH/mc is the cyclotron frequency and I =j— ij/ Dr. 

We turn now to the interaction correctionsEl-to the 
conductivity. Our analysis follows Choi et aLcil Note 
the interaction MR is mainly due to the magnetic 
field dependence of the pre-factor of these corrections, 
f{Si,Sh,Q,H) in Eq. (|l), rather than the explicit de- 
pendence upon H expressed in Eqs.p^^^||_|||). 



The diffuson channel correction i: 



(4) 



5oln(^) 



for ^ > 1 



.2-f^5o[*(i + ^)-*(5)]-&2W, for^<l 
where 

90 = 4-3-^ln(l + -j, 



2 + F 
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dO 
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ln(l + |)-4, 
1 H sm I - 



(5) 
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h 



/ dfi—^ In 

Jo 

g*^iBH/T, 



and 



, . _ r ln(Vl-3) for /i> 1, 
9^^ ' \ 0.084/i2 for/i<L 

\l/ is the Digamma function, k is the 2D inverse screening 
length, i/(f2) is the density of states at an energy hft, g* 
is the 2DHG effective g-factor, and is the Bohr mag- 
neton. In the diffuson channel the interaction correction 
depends on the magnetic field through the Zeeman en- 
ergy only. The diffusive correction, Sa^^f, is hence inde- 
pendent of H as long as /i <C 1. Unlike the cooperon 
correction that vanishes for lOcT > 1/g, the diffuson cor- 
rection persists to UcT ~ 1. 
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T„(S,)=4.9+.2K , T„(S^)=5.8+.2K , T„(Q)=5.1 + .2K , T„(p )=4.3+.2K 
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FIG. 3: The various scattering rates expressed as resistances 
(left axis) and the zero field longitudinal resistance (right axis) 
vs. T. Solid lines depict best fit to Arrhenius dependence with 
the characteristic temperatures listed at the top of the figure. 
Inset - same data in semi log plot. 



The cooperon channel correction to the conductivity is 
given byO 



where 



tdt 



2eDH\ 
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sinh^(t) L sinh(a;i) 

ln(a;) for a; 3> 1, 
0.30x2 fQj. 2; < 1, 



and 



g:\T,H) = 



9o 



i+ln(Li3^) for2|0|L<l, 



9o 



for 



2eDH 
ttTc 



> 1. 



The correction depends on H and vanishes for H > 
27r_fftr, or equivalently, \/2^cTg/27r > 1. The Maki- 
Thompson correction is ignored. 



C. Data analysis 

For the high mobility sample, the WL and WAL cor- 
rections are negligible above 1 Gauss. The analysis may 
thus be confined to interband scattering and interaction 
in the diffuson channel [Eol (|)EI]. The full MR is de- 
picted in the inset to Fig. g. The parabolic background 
is attributed to Coulomb interaction and analyzed be- 
low. Subtracting this background we obtain the data 
presented in Fig. 0. Fitting these data to Eqs. (|2|), one 




FIG. 4: (a) High mobility sample. ASi ^ Si - 2.16Q, 
ASh ^ Sh- 2.16-1(3 vs. T for p = 2 X lO^cm^^. (b) Hole- 
hole interaction correction to the conductance in the diffuson 
channel, Saf^t- 



obtains Pxx{H — > oo), L, and as a function of tem- 
perature and density. Note the excellent agreement with 
the predicted shape, Eqs. (||). Below ~ 0.6 K, the resis- 
tance is practically independent of T while for tempera- 
tures above ~ 2 K, the Lorcntzian is hardly visible. The 
suppression of the classical two band magnetoresistance 
results from interband scattering. At low temperatures 
this scattering is mainly elastic. As the temperature is 
increased, inelastic scattering commences, the drift ve- 
locities of carriers in the two bands gradually approach 
each other, and the magnetoresistance is consequently 
diminished. The extracted functions. Si, Sh, and Q, to- 
gether with Pxx{H — 0), are shown vs. T in Fig. || for the 
same total density as in Fig. ||. At low temperatures, all 
these quantities saturate to some residual values S'((0), 
5/1(0), (5(0) which we attribute to inter and intraband 
elastic scattering. As the temperature is increased, in- 
elastic scattering commences and these quantities grow. 

The remarkable and central result emerging from the 
high mobility data is the observation that the inelastic 
scattering rates follow the same temperature dependence 
as Pxx{H = 0), namely, S^{T) = S^{0)+as, exp{-To/T), 
Q{T) — Q(0) -I- aq exp{—To/T), where 05., ag are con- 
stants. The characteristic temperature, Tq, is similar to 
all these quantities, including pxx (see inset to Fig. ||). 
For Si we find Tq = 4.9 ±0.2 K, for Sh we find Tq = 5.8 ± 
0.2 K, for Q we find Tq = 5.1 ±0.1 K, and for pxx we find 
To = 4.3 ± 0.1 K. We find experimentally [see Fig. |(a)] 
that for T < 1.6 K, Si{T) ~ Si{0) + 0A6-^[Q{T) ~ Q{0)]; 
Sh{T) ~ Sh{Q) + 0.46[g(r) - Q(0)], thus yielding for 
the density of Figs. || and ||, a = 0.46. Since the pre- 
factors multiplying Q are reciprocal, the resistances 5/(0) 
and Sh{0) are identified as the diagonal resistances, per- 
taining to intraband scattering. The Arrhenius tem- 
perature dependence is hence traced to inter-band scat- 
tering alone, rendering the intra-band scattering practi- 



FIG. 5: Solid lines - heavy and light particle-hole excitation 
continua as a function of momentum scaled to the heavy hole 
Fermi wave vector. Shaded area corresponds to the range 
where drag-like interband scattering is possible at very low 
T. Dashed lines -optical (op) and acoustic (ap) plasmon 
dispersions. Inset - The measured bands dispersion rela- 
tions. The energy Ef corresponds to the hole Fermi energy 
at ptotai = 4.25 X 10" cm-^ 



cally temperature independent. We attribute this unex- 
pected scattering rate to the acoustic plasmon mediated 
Coulomb interaction discussed below. The results are 
identical to those presented in Ref. 

We turn now to analyze the interaction correction to 
the conductivity. The Coulomb interaction is character- 
ized by F, the angular average of the statically screened 
Coulomb interaction, given by Eq. (|). With Si,Sh,Q 
known from the previous fit, we fit the parabolic back- 
ground to Eq. (0) to extract Saint- We find F = 0.5 ±.05 
for all T except the lowest ones. This value should be 
compared with the theoretical value, F — 0.81. Fig. ||d 
compares the experimental Saint with the theoretical one 
as a function of T. For temperatures above 0.8 K, the 
agreement between theory and experiment is excellent. 
At lower temperatures the experimental interaction cor- 
rection is smaller than that predicted by theory. The in- 
teraction coefficient, go{F = 0.47) [Eq. (|)], is 0.65±0.03, 
compared with go{F = 0.81) = 0.46, predicted by theory. 
It should be noted that the theoretical value is calculated 
for two identical parabolic bands with m — 0.3mo. The 
actual band structure is more complex and in fact, non- 
parabolic. The fact that F is independent of temperature 
above 0.8K supports the attribution of the interaction 
contribution to the diffuson channel. We emphasize that 
the saturation of the interaction correction below 0.8K 
is real and does not result from carrier heating. The lat- 
ter, deduced from the SdH data, proves efficient cooling 
down to 0.2K. 

In section ^ we discuss in detail possible reasons for 
the Arrhenius temperature dependence of Si , Sh and Q 
which in turn leads to a similar temperature depen- 
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FIG. 6: Low mobility sample. Longitudinal resistance vs. 
temperature at different carrier densities. In the metallic 
regime the data are fitted to an Arrhenius function. The 
characteristic temperature, To, is indicated in the figure. 



FIG. 7: Low mobility sample. Longitudinal resistance vs. 
temperature for different parallel magnetic fields. PMR and 
suppression of the metallic characteristics by the parallel field 
are evident. 



dence of pxx- These scattering rates (expressed as re- 
sistances) crucially depend on the bands' dispersions re- 
lations, i?i(k), and their resulting excitation spectra. To 
extract the bands dispersions depicted in the inset to 
fig. ^ we approximateEflLJ the light band by a parabolic 
relation with a mass mi = 0.28too (mo is the bare elec- 
tron mass) which was found from SdH temperature de- 
pendence. The variation of p;, ph with ptotai, depicted 
in fig. |l], is then used to calculate the ratio between the 
two bands compressibilities. Neglecting band warping as 
well as differences between density of states and com- 
pressibility, we use the ratio of the two compressibilities 
to extract the dispersion of the heavy band. This dis- 
persion then allows, within the random phase approxi- 
mation, the calculation of the excitation spectrum of the 
system. The spectrum is composed of two particle-hole 
continua, one for each band, and two plasmon branches. 
Both are shown in fig. ^ for zero temperature. 

In summary, we find that the MR of high mobility 
2DHG is governed by classical, two band PMR at weak 
fields and NMR (parabolic background) due to Coulomb 
interaction in the diffuson channel at stronger fields. The 
resistance increase with temperature ("metallic phase") 
follows the inelastic interband scattering. Careful anal- 
ysis of the interaction contribution to the conductance 
reveals unexpected saturation of the interaction below 
0.8K. 



III. LOW MOBILITY 2DHG SAMPLE 

The transition to an insulating behavior in the high 
mobility sample is taking place at a very low carrier con- 
centration where the sample becomes inhomogeneous. 
We have therefore preferred to study the transition 
regime and the insulating phase in a lower mobility sam- 
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FIG. 8: Low mobility sample. Normalized magnetoresistance 
at different carrier densities. Data points are marked with 
crosses while lines correspond to the theory discussed in the 
text. Note the crossover from PMR to NMR at a density just 
above the MIT. 



pie, ^ ~ 20,000 cmV(Vs) for p = 4 x 10" cm^^ ^t 
400 mK, where the transition occurs at a moderate den- 
sity. 

Following the high mobility data analysis, we study 
the magnetoresistance of the low mobility sample and 
relate it to the temperature dependence of pxx at H = 
0. Since the quantum corrections to the conductivity 
are order unity in quantum conductance units, the WL, 
WAL, and interaction correction are more pronounced 
here compared with the high mobility sample. 

In Fig. H we present normalized MR measurements 
at 400 mK for various densities. For high densities, 
p > 1.5 X 10^^ cm^^, the weak field MR is positive 
and dpxx{T,H = 0)/dT > 0. For 1.5 x 10" cm'^ > 
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FIG. 9: Low mobility sample. Magnetoresistance curves 
at different temperatures. Inset - magnified view at T = 
0.4, 0.9, 1.3, 1.5, 2K. Note the low T, weak field cusp charac- 
teristic to WAL. 



FIG. 10: Low mobility sample. Subtraction of the quantum 
correction from the longitudinal resistance according to Eq. ^ 
yields the classical Lorentzian PMR. 



p > Pc = 1.36 X 10^^ cm^^ the MR is negative while 
dp/dT is still metallic. For p < pc both the MR and 
dpxx {T, H = 0)/ dT are negative. From the two SdH fre- 
quencies one may extract the carrier densities, pi andp^i, 
for the lighter and heavier bands respectively, as depicted 
in Fig. ||. Below ptotai ~ 2 x 10^^ cm~^ the two bands 
are no longer distinguishable. At higher densities, the 
expected Lorentzian PMR is found. The MR curves re- 
mind Bergmann's dataE3 on Mg thin films covered with 
a fraction of atomic gold layer. In those experiments a 
crossover from WAL to WL was found as function of the 
spin-orbit scattering strength tuned by the amount of 
gold. In our experiment, the role of gold is played by the 
energy gap between the spin-orbit split bands combined 
with elastic interband scattering. 

We focus first on the high density regime. A set 
of resistance curves vs. magnetic field, for ptotai — 
4 X 10"" cm-2, pi = 1.7 X 10"" cm'^, p^ = 2.3 x 
10"^^ cm~^ and different temperatures is depicted in 
Fig. ^. The PMR is dominated by classical MR but near 
the origin, as evident from the inset to Fig. the exper- 
imental points deviate below the Lorentzian curve and 
form the well known WAL cusp. Eq. (]l|) contains six in- 
dependent variables; Pl{H — > oo), L, W , F, and Tso- 
To find them we take advantage of the different magnetic 
field scales characteristic to WL and WAL compared with 
the classical MR. Using H > lOHtr ^ 500 Gauss data 
and ignoring WL or WAL corrections we determine 4 
parameters, pl{H oo), L, W, and F. We then re- 
sort to small magnetic fields and find the remaining two 
parameters, and Tjo. 

Though the procedure for determining the four param- 
eters is similar to the one used in the high mobility case, 
the larger quantum corrections lead to uncertainty in F 
and consequently, in Sh- We have therefore adopted the 
following strategy. We first fix F to its theoretical value. 
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FIG. 11: Low mobility sample. ASi and AQ (left axis) and 
the zero field longitudinal resistance (right axis) vs. T. Solid 
line depicts best fit to Arrhenius dependence upon tempera- 
ture. 



F — 0.8, and calculate Pl{H oo), L and W which 
in turn yield Si, Sh, and Q. However, the resulting Sh, 
and Q display a slight nonmonotoneity upon tempera- 
ture below 0.8 K. Then, for T < 0.8 K we tune F to give 
a monotonous increase of Si, Sh, and Q with T. The 
needed change in F is less than 10%. Different values of 
F, {F = 1,0.5) give similar results with different Sh- For 
F = 1, the overall change in F is minimal (5%). Fig. |l^ 
depicts the MR data for F = 1 after subtraction of the 
quantum corrections. The agreement with the expected 
Lorentzian is very good. In Fig. |l|. Si -Si{T^ 0.4 K), 
Q - Q{T = 0.4 K), and p^xiH = 0) are depicted vs. T 
for the same total density as in Fig. g and F = 1. Unlike 
the high mobility case, the uncertainty in F and hence, 
Sh, might be substantial. Consequently we limit our- 
selves to a qualitative comparison between Si, Q , and 
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FIG. 12: Low mobility sample. Left axis: Hole-hole inter- 
action correction to the conductivity in the diffuson channel, 
Sa^nf Right axis: Angular average of the statically screened 
Coulomb interaction, F. Theoretical value is F = 0.9. 



Pxx{H = 0). Examining Fig. ITT] one finds a qualitative 
agreement between the temperature dependence of tlie 
three resistances, in agreement with the high mobihty 
data. 

Fig. |l2| compares the theoretical Saint (for F = 1), 
with the measured one as a function of T. For temper- 
atures above 0.8 K, the agreement between theory and 
experiment is satisfactory, but for temperatures befow 
0.8 K the interaction correction is again smaller than the 
predicted value. We emphasize that theory was done 
for two parabolic bands with m = O.Smg. The pro- 
nounced non-parabolicity of the heavy hole band may 
affect the interaction contribution significantly. Unlike 
the high mobility case we do not observe saturation at 
low temperatures. From the mass analysis of the SdH 
oscillations, depicted in Fig. ^ one finds mi — 0.28mo, 
and no carrier heating above 0.3 K. 

Once Pl{H — > oo), L, W, and F are determined one 
can fit the low magnetic field data to <j{H) — a{H = 0) = 
Saa-L^H) + SamtiH) - Sa^vLifi) - 5aint{^) and extract 
Tip and Tso- Since the Lorentzian width, is larger 
than Htr, and since the ratio of the thermal and de- 
phasing lengths make the weak field interaction contri- 
bution ^ 5_g^ times smaller than the WL or WAL one 
[see Eqs. (||), <@], the MR for H < Htr is dominated 
by WL or WAL. The uncertainties in F are hence re- 
flected in less than 10% uncertainty in and Tso- Fig. |lj 
presents a{H) - a{H = 0) - [aci{H) - 0-^(0)] for ptotai = 
2.8 X 10^^ cm~^ where crci{H) is the classical Drude part 
of the conductivity. The data agree well with theory, 
Eq. (^). The resulting and t^o, forp — 2.8 x 10^^ cm^^ 
are depicted in Figs, p^and |l6|. The error bars reflect the 
uncertainty in F. The phase breaking rate is linear with 
temperature, as expected by theory (Fig. p^ ). The pref- 
actor, however, is five times larger than the theoretical 
prediction 1/r^ = T/{hg) ■ hi{g/2). These r^p values are 



FIG. 13; Shubnikov-de Haas mass analysis. Best fit yields 
mi — 0.28mo where mo is the bare electron mass. 
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FIG. 14: Low mobility sample. a{H)~a(0)-{aci{H)-aci{0)) 
vs. T for ptotai = 2.8 x lO^^cm"^. aci{H) is the classical 
Drude conductivity. Theoretical curves correspond to eq. 
Note the PMR in the whole WAL field range. 



consistent with those found in n-type SiliconEj and p-type 
GaAsEJ systems. The spin-orbit scattering time, Tgo, is 
as expected practically independent of temperature. The 
extracted times satisfy > Tso and comparable to 
T. Consequently, the WAL PMR should persist in the 
whole range, < H < Htr- The band splitting, eg, may 
be roughly estimated with the help of the Dyakonov-Perel 
formula and compared with the value extracted from the 
dispersion relation (inset to Fig. ^ deduced from SdH 
oscillations. Close to the MIT the Dyakonov-Perel spin 
precession depends on the relaxation time, r, extracted 
from the decay of the SdH oscillations once the mass has 
been determined. Fig. |l3 depicts the SdH mass analysis 
for p = 2.8 X 10^^ cm . The envelope analysis yields 
T ~ 10~^^ sec. Substituting r in Dyakonov-Perel ex- 
pression one finds eg ~ 0.5 mev, which agrees with the 
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FIG. 15: Low mobility sample. WAL analysis, (a) dephasing 
rate vs. T at different densities, (b) spin-orbit scattering time 
vs. T at different densities. 

SdH value. This result strongly supports the association 
of the WAL with scattering between the two spin-orbit 
split bands and establishes WAL as a powerful tool for 
measuring eg. 

We turn now to analyze the data corresponding to den- 
sities below p = 2 X 10"'^^ cm~^ where the Lorentzian 
classical MR and the two SdH frequencies can npt be 
resolved. Then, Eq. (|l|) takes the single band formcll, 

PxxiH) = po- Po[l- (wct)^] {Sa^vL + Saint) , (6) 

where po is the classical longitudinal resistance. There 
are 4 parameters to be determined, po, F, and Tso- 
Again, we separate the fitting procedure to 3> Htr 
and H ^ Htr- Albeit, this time Htr > 0.1 T and conse- 
quently the WAL or WL corrections are important in the 
whole magnetic field range, < H < 0.7 T. The four pa- 
rameters are found iteratively. In the weak field regime 
we guess po and F and fit a{H) — a{H — 0) to theory 
with and Tso as two fitting parameters. The resulting 
and Tso are plugged into Eq. and the high field 
data are fitted with po and F as two adjustable param- 
eters. This procedure converges after ~ 10 iterations to 
a self consistent solution. The results of such an analy- 
sis for p = 1.5,2 X 10^^ cm^^ are presented in Figs. |l^ 
and |l6|. Again, Tso is practically independent of tempera- 
ture. Fig. ^ depicts the quantum correction to the mag- 
netoconductivity (j{H) — a{H = 0) — [aci{H) — (Jci{0)\ 
for ptotai = 2 X 10^"'^ cm~^. The data agree well with 
Eq. (|). Since H^o < Htr the MR in Fig. |l| changes 
sign from positive to negative as the field increases. The 
weak field PMR results from WAL {t^ > Tso) and one 
may use the Dyakonov-Perel formula to estimate eg . For 
p = 2 X 10^^ cm^^ and p = 1.5 x 10^^ cm^^ we find 
eg ~ 0.3 mev ~ 0.2Ep and eg ~ 0.2 mev ~ O.IGE^, 
respectively. For p > 2.8 x 10^^ cm~^, the classical 
PMR as well as the two frequencies of the SdH oscilla- 
tions are visible implying eg > Ti/t. For densities below 



FIG. 16: Low mobility sample. Tlp/tso vs. total carrier den- 
sity (left axis), t/tso vs. total carrier density (right axis). 
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FIG. 17: (j{H) - cr(O) - {(Ja[H) - a^iiQ)) vs. T iovptotai = 
2x10^^ cm~^. Theoretical curves correspond to Eq. (0). Note 
the crossover from PMR to NMR resulting from the fact that 

T < Tso. 



2.8 X 10"'^"'^ cm~^, eg < Ti/r, and the two bands can be 
resolved solely by WAL. These results demonstrate that 
the Lorentzian MR as well as the two SdH frequencies 
disappear when the band spacing becomes smaller than 
the level broadening by disorder. The quantum interfer- 
ence, on the other hand, indicates the existence of two 
bands down to lower densities. 

The extracted Coulomb interaction parameter is ~ 
1 ± 0.2, practically independent of temperature. This 
value is in good agreement with the theoretical predic- 
tion, F — 0.9. The resulting po increases with tempera- 
ture. The origin for this temperature dependence is not 
clear. It might result from inelastic inter-band scattering 
or other effects such as percolation or temperature de- 
pendent screening. Those mechanisms are discussed in 
the next section. 
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FIG. 18: Low mobility sample. a{H) - a(0) - [aci{H) - 
aci{0)] vs. T ior ptotai ~ 1.36x10^^ cm~^. Theoretical curves 
correspond to the standard WL formula [Eq. (^) with Tso —* 
oo]. 

For p < Pc = 1.36 X 10^^ cm^^, WL and interaction 
dominate and lead to NMR. At these densities, g < 2.5, 
the quantum corrections are on the order of the classical 
resistance, and the expansion in 1/g is only approximate. 
There are three parameters to be determined: F, and 
T. Fitting cr{H) — a{H = 0) to theory indicates that WL 
dominates over interaction. Different values of F modify 
the extracted by less than 10%. One may therefore 
set F = 1, and find and r. Fig. O depicts the data 
and a best fit for ptotai — 1-35 x 10^^^ cm^^. Fig. |l^ 
displays the extracted dephasing rate for the three lowest 
densities. Note that extrapolates to a finite value 

for T 0. Similar saturation of has been reported by 
Brunthaler et alM in their Si-MOSFET weak localization 
measurements. We emphasize though that we have used 
WL theory at the border of its validity, Sg < g. The 
extracted values might hence be misleading. 



IV. DISCUSSION OF DATA ANALYSIS 

At high densities, p > 2 x 10^^ cm^^, we find that 
the PMR and temperature dependent resistance are fully 
consistent with scattering between the two spin-orbit 
split bands, WAL, and interaction. The resistance in- 
crease with temperature is mainly a consequence of in- 
elastic interband scattering. For intermediate densi- 
ties, 1.36 X 10" cm-2 < p < 2 X 10" cm^^ WAL 
indicates two non-degenerate bands. For example, for 
p — 2x 10^^ cm~^ we find eg/ Ep ~ 0.2, implying roughly 
20% difference in densities and resistances of the two 
bands. For densities smaller than the critical density, 
Pc = 1.36 X 10^^ cm~^, the two bands broadened by dis- 
order merge to form a single band. The classical PMR as 
well as the WAL vanish and leave a conventional, single 



FIG. 19: WL analysis of low mobility sample. Dephasing rate 
vs. T at different carrier densities. Note the deviation of t^^ 
from linearity in T as discussed in the text. 



band WL NMR. Similar MR is alsOrioprid in other ex- 
periments on p-type GaAs and SiGec3G3 samples. The 
zero field WL crosses over to strong localization as the 
density is further reduced. The MIT we observe reflects 
in our opinion effective cancellation of the temperature 
dependence of the resistance due to two competing pro- 
cesses. On one hand the reminiscent band splitting leads 
to inelastic interband scattering and WAL, both charac- 
terized by dpxx/dT > 0. On the other hand, interaction 
and WL lead to dp^x/dT < 0. As the density is re- 
duced, the splitting between the two bands shrinks, their 
broadening due to disorder is enhanced and they prac- 
tically merge to form a single, doubly degenerate band. 
In that process, the effect of inelastic scattering on the 
resistance is gradually diminished, the WAL turns into 
WL (r^ becomes comparable to r^o), and together with 
the Coulomb interaction, the sample displays an insulat- 
ing behavior. The so called MIT is hence just a smooth 
crossover from two to one band physics. Note that on 
the metallic side, |ii^iaij the transition, our MR data, as 
well as other data,c2lE3 agree well with WL. Nonetheless, 
the temperature dependence is opposite to that expected 
from WL. There should, hence, be some scattering mech- 
anism with opposite T dependence compared with WL 
and very weak dependence upon magnetic fields. Extrap- 
olating our knowledge from the metallic regime we con- 
jecture that this mechanism is inelastic interband scat- 
tering. 

In other, .experiments on higher mobility 2DHG 
samplesBll3o, the critical density, Pc, turned out to be 
lower than that reported here. In Ref. ^ pc = 1-25 x 
10^*^ cm~^, and in Refs. p^ , ^ , pc is around 5 x 10^° cm^^. 
We believe that in those experiments the bands are split 
down to densities lower than in our experiment, due to 
the following reasons: (a) For the mobilities of Refs. pyiC|, 
the bands broadening by disorder, and hence the bands 
mixing, is small, (b) The quantum well asymmetry 
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and hence the lack of inversion symmetry are more pro- 
nounced in those samples. Consequently, the band split- 
ting due to spin-orbit interaction is considerably larger 
than in our samples even with the lower densities taken 
into account. The lack of inversion symmetry is partic- 
ularly pronounced in the p-type invertednSemiconductor 
insulator semiconductor (ISIS) structureEa used in Ref. ||. 

Quantitative estimates of eg require detailed band 
structure calculationsLJ. WinklcrEZi studied band split- 
ting in 2DHG due to lack of inversion symmetry and was 
able to show that SdH analysis underestimates the differ- 
ence in carrier partition between the two bands at = 0. 
For a typical confining potential at p = 2 x 10^° cm^'^, 
Winkler found a 10% difference in hole densities and 20% 
difference in their masses. The splitting, eg, turned out 
to be about 0.5 K. These results confirm the role of two 
bands in the transport properties of low density 2DHG. 

We turn now to discuss some other aspects otthe sam- 
ples that might be relevant for the MIT. Meircll pointed 
out that the 2DHG is probably inhomogeneous near the 
MIT and analyzed our data in the context of classical 
percolation. He was able to show that the H — low 
temperature longitudinal resistance can be fitted with 
p^^ = a{p - p*)-^, with p* = 0.9 ± .05 x 10" and 
7 = 1.3 ± 0.1. For classical percolation, the critical 
exponent 7 should be 4/3 ~ 1.33. Identical functional 
dependence was found for the data of Hanein et al£3 In- 
dications of inhomogcneity are found in thermodynamic 
measurements as well. Dulta-|and Jiang measured the 
compressibility, k, of 2DIIGE2I, and found that dn/dp 
changes sign at the critical density. We also measured 
the sapie quantity in a different way and found similar 
resultfM. Si and Varma argue that the inverse compress- 
ibility should-vanish as the MIT is approached from the 
metallic sidcEJ. The argument is that near the transition 
disorder and Coulomb interaction create inhomogeneous 
regions (puddles) which behave like weakly coupled quan- 
tum dots. The breakdown to puddles leads-to Coulomb 
blockades and incompressibility. Hani et alEd used a sin- 
gle electron transistor to measure the local compressibil- 
ity of a 2DHG and found that below the critical MIT den- 
sity, the sample is indeed no longer homogenous. At the 
moment, the importance of inhomogcneity for the MIT 
is unclear. Inhomogcneity should have an interesting ef- 
fect on the phase breaking length, l^. When puddles are 
formed, the diffusion within each puddle is much faster 
than that between puddles. In that case one may imag- 
ine a situation where for a limited temperature range 
seemingly saturates to the puddle size. Fig. ^ depicts 
= D/l'^ (where is extracted from WL and D is the 
diffusion constant extracted from the global resistance) 
as a function of temperature. Indeed naively extrap- 
olates to a non zero value for T ^ 0). We do not have 
though any evidence that links the deviation of t^^ from 
linear T dependence with puddle formation. 

Another source of a metallic behavior of 2D systems 
has been proposed—by Gold and Dolgopolov, and Das 
Sarma and Hwangcil, and was very recently revisited by 



Zala et allB These authors consider impurity screening 
by the 2DEG. The main dependence of the dielectric 
function upon temperature results from 2kp back scat- 
tering of the electrons. The calculation suggests that 
Aa/(T — f{T/Tp) where / is some scaling function. The 
disorder potential in Si-MOSFET is short ranged and 
characterized by significant 2kp wave- vector components. 
In modulation doped semiconductors, however, the sepa- 
ration of the dopant layer from the 2D gas leads to a long 
range impurity potential with wave vectors considerably 
smaller than kp. Nonetheless, Hamilton et alEB took the 
function f{T/Tp) as universal function and succeeded to 
collapse their GaAs hole data (various densities), to a sin- 
gle curve which is linear at small T/Tp, and saturates to a 
constant value at higher values. We tried a similar analy- 
sis but. could not find such a universal scaling. Hamilton 
et alBB argue that slightly above the critical density, the 
asymmetry of the confining potential is irrelevant for the 
resistance change with temperature. By fixing g = kpl, 
they claim to eliminate the influence of the confining po- 
tential shape and conclude that the metallic behavior is 
due to temperature dependent screening rather than two 
bapds physics. However, in a recent paper, Papadakis et 
al.t3 show that in less symmetric quantum wells, where 
subband splitting is large, the resistance increase with 
temperature is large and uncorrelated to the low T re- 
sistance. While temperature dependent 2kp screening 
might be a relevant factor, we find the two band physics 
essential for understanding our results. 



V. THEORETICAL ASPECTS OF 2DHG 
MAGNETOTRANSPORT 

A. Arrhenius temperature dependence of the 
hole-hole interband scattering rate 

1. Qualitative discussion 

In this section we consider the hole-hole scattering con- 
tribution to the resistivity in a two band system. Since 
the source of band splitting in the system is spin-orbit 
coupling, two momenta states in different bands gener- 
ally have overlapping spin wave functions. A finite in- 
terband scattering at vanishingly low temperatures (see 
fig. ^) results from impurity scattering. At higher tem- 
peratures, hole-hole scattering commences and due to the 
different mobilities in the two bands, increases the resis- 
tance. Similarly to the H = Coulomb drag between 
two layers, one naively expects a Coulomb scattering 
contribution to the resistance. Experimentally we rather 
find an unexpected Arrhenius dependence upon tempera- 
ture. Similar T dependence is found in many experiments 
probing the metallic phase in 2D. 

We turn now to analyze the temperature dependence 
of hole-hole inter-band scattering and show that its en- 
hancement by acoustic plasmons leads to an Arrhenius 
law for the resistivity at temperatures which are not too 
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FIG. 20: The {q,u}) plane for small q and u. The regions in 
which interband and intraband particle-hole excitations are 
possible are marked in the figure. The filled regions mark the 
intraband excitations in the light band {££) and the interband 
excitations from the heavy band to the light band {h£). The 
scattering process discussed in the text [Eq. (^)] is limited to 
q and uj in the overlap of both regions. 



low. 



As discussed in Section V A 2 below, the dominant con- 



tribution leading to an Arrhenius T dependence arises 
from scattering of a light hole off another hole that 
changes its band in the process (light to heavy or vice 
versa) . The contribution of this process to Se is given by 



1 h 1 



47r2 e2 Tpj J {2n 



smh\huj/2T) 



(7) 



X \ysc{q, ImXw (q, ^) irnxwl?, ^) , 

with similar expressions for the contributions to Sh and 
Q. In Eq. V^c(q, i^) is the screened Coulomb interac- 
tion. The X functions are response functions generalized 
to a two band system. They are defined explicitly in 
Eqs. (p^) and (p^. Their imaginary parts, appearing in 
Eq. (Q)7are non vanishing whenever it is possible to ex- 
cite a particle-hole pair with momentum ?iq and energy 
Tiuj by exciting a particle from band i to band j (with i, j 
being l^ l or /i, I). 

Fig. ^ depicts regions in the (g, lo) plane for which 
particle-hole excitations of different types exist (for small 
q and lo). Intraband excitations are possible only to the 
right of a line whose slope at the origin corresponds to the 
Fermi velocity in that band (electron- hole continuum). 
Interband excitations are forbidden at very small q and w. 



At w = the minimal wavevector transfer is k 



while at g = 0, the minimal energy for excitation from the 
heavy to the light band corresponds to the band splitting, 
Eg, between the bands at kp . The scattering process 
described by Eq. (|^) is possible at q and lo for which 
both Imxu and Imx^^ are nonvanishing. The various 
scattering processes are illustrated in fig. pTl. 




FIG. 21: The arrows represent the different particle-hole 
excitations. This figure is given for illustration only — the 
real situation is more complex as momentum space is two- 
dimensional. 




FIG. 22: The (g, lu) plane, for small q and lu. The optical (OP) 
and acoustic (AP) plasmon branches are shown in dashed 
lines. The regions in which interband and intraband particle- 
hole excitations are possible are marked as in fig. pol The 
scattering process discussed in the text [eq. (^)] is strongly 
enhanced along the acoustic plasmon line. As seen in the 
figure, there is a minimal frequency, coo, for this process. 



Due to the existence of two bands, the system is char- 
acterized by two plasmon branches. The optical branch is 
essentially governed by the motion of light holes. Its dis- 
persion relation is similar to that of plasmons in a single- 
band 2D system, namely, lo oc ^/q. The acoustic branch is 
characterized by a linear dispersion, lo (x q, correspond- 
ing to the motion of the heavy holes whose mutual inter- 
action is screened by the light holes. The acoustic mode 
is hence similar to a sound wave in a metal, namely, ion 
plasma oscillations screened by electrons. The ions are 
played here by the heavy holes while the electrons are 
replaced by light holes. The acoustic plasmon velocity is 
approximately given by Vp = -\/to^/2to£ vi!^\ certifying, 

(h) ^ ^ {i) 
V}. ' < Vp < V},'. 

The Coulomb scattering depends on the screened 
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Coulomb potential, Vsc{q,uj) = Vb{q)/£{q,uj), where 
e{q,oj) is the dielectric function of the system. A large 
contribution is therefore expected for q and w close to 
the plasmon branches where e{q,uj) — 0. The Arrhe- 
nius temperature dependence of the interband scattering 
rates results from plasmon enhanced scattering. A simi- 
lar effect was previously considered for Coulomb drag in 
double- layer systemsEi . 

In fig. ^ we plot the {q, u) plane, this time with the 
plasmon dispersion lines added. Plasmon enhanced scat- 
tering occurs for the process described above only above 
a threshold frequency luq- For temperatures T hioo, the 
main contribution comes from frequencies just above the 
threshold, ujq < uj < ujQ + T/h. The contribution of these 
scattering events is suppressed by a factor c>c e~^^l'^ due 
to phase space considerations but the plasmon enhance- 
ment, compensates for that reduction. The resulting 
temperature dependence is Arrhenius, cx ^^^^oIt ^ 

We turn now to a detailed calculation of the intra 



and interband scattering rates resulting from plasmon 
enhanced hole-hole scattering. 



2. Detailed calculations - Boltzmann theory for two bands 

We start by solving Boltzmann equation for the hole- 
hole interaction contribution to the inter and intraband 
scattering rates. 

For a space and time independent distribution func- 
tion, /(*-'(k) {i — l,h), the Boltzmann equation in the 
presence of an electric field E takes the form 



eE df^^_(df^^' 



dt 



(8) 



The collision integral for hole-hole scattering is given by 



h-h ji'j' k2k'jk^ 



1 - f^'^ 



;'ik2^kik2 J 



(9) 

where Wl-'fC! is the scattering rate for two holes in states ki,i and k2,j to scatter to and k2, j' . Hole-hole 



kik2^kjk2 



scattering satisfies both energy conservation, W^k{k*2*-ik' k' ~ ^ki^2^k' k' ^ (^fci 
conservation, ^'^kik2^k'k' '^ki+k2,k'j^+k^- Using the golden rule we have 



4'^ - ^'^ 



and momentum 



S^k,-q,k2+q = Y|M,Hki, ki - q)|2 |Af,,,(k2, k2 + q)P \VUq,^) 



(10) 



Here, Mij(k, k') are the matrix elements between the spin wavef unctions, derived in Section VC, below. 



For convenience wc take = Kk/rm, i.e., neglect non-parabolicity and non-isotropy of the two bands. Linearizing 



the collision integral for hole-hole scattering, we have 



/a/W(ki) 

^ dt 



h-h 



fo[^^]-fo{e-'-u.]\S(e-'-e 



jl'j' ksq 



hduj 



-oo Smh\hL0/2T) klk2-ki-q,k2+q 



f ii') 



.p') 
k2+q 



(11) 



'fel 



1kl-q| 



— Hlu 



(i) IJ'} 



1k2+q| 



We define the intra and interband response functions, 

/o (4*^) - fo (4k+q|) 



Xij 



iq,i^) = -J2- 



k "=fe 



Xi) Xj) 



«|k+q| +^^ + ^^ 



Afy(k,k + q) 



(12) 



In addition, we define response functions with momentum (or current) vertices by adding an appropriate function of 
k to the definition 



s(k) 



Xj) 
1k+q| 



|M„(k,k + q)|^ 



(13) 



For example, the current-current intraband response function is written in this notation as ^^^(q, w), where j 
(k -|- q/2)/mi. In this example, the function g(k) depends on q as well. 
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Solving the linearized Boltzman equation we derive the contribution of hole-hole scattering to the resistivity matrix, 
1 h I r°° Muj 



h-h 
Pik 



EE 



47r2 e2 TpiPk 



smh\huj/2T) 



\Vsc{q.io)\' 



I (_q^ -c^)Inix,,' (q, io)6,k + Irax'^-, (-q, -c^)InixJj+'== (q, uj)6j,k 

- ImXii' (-q, -'^)ImXjj' (q, - Inix,,/ (-q, -uj)lraxfy (q, L^)6jk j 



r 



(14) 



The dominant contribution to Eq. is given by 

Eq. (Q). In this process a hole changes its band by scat- 
tering off a hole in the light band (i = j = j' = i, i' = h 
OT i = i' = j' = £, j = h OT i = i' = j = £, j' — hX Other 
contributions are readily extracted from Eq. (|lj). For 
example, for the ££ component of the resistivity matrix 
there are four additional contributions: 

• Interband Coulomb drag (i = i' — £, j = j' = h). 
This is the only term remaining in cases of van- 
ishing matrix elements between the spin wavefunc- 
tions of the two bands, e.g., in the absence of spin- 
orbit coupling: 



1 



h duj 



h 1 

X |V;c(<?,w)pImx«(g,w)ImXh/i(g,w 



(15) 



.(2) 



Band exchange (z = j' ~ £, i' = j = h): 



1 h 1 

"T2 12 rr„2 



47r2 e2 Tpj J {2t:)^J_^ smh\huj/2T) 



X <^ ImX;j?(q,w)ImxM(<7,t^) - Imx^|(q,w) 



Vsciq,to)\' 



(16) 



• Two holes scatter from one band to the other (i = 

j = £, i' = j'^hy. 



.(3)- 



1 h 1 r d^q 
A^'^Tpj J {2ttY 



hduj 



smh\huj/2T) 



l^sc(g,a;)|' 



: {imxSq, c.)Imx..(<Z, c.) -f Imx,^;^(q, .)Imx^/-(q, .)} . 

(17) 

• One particle changes its band by scattering off a 
heavy hole {i = £, j — i' — j' — h): 



.(4) 



I h I r d^q hduj 

4^ ^ I>f y (2^ sinh2(?ico/2T) 



X |v;c(g, t^)P Imx^^ (q, ^) imxhhiq, t^) • 

The analogous expression for p^h is derived directly 
from pii by interchanging band indices, £ h. Using mo- 
mentum conservation, we find that the general structure 



of the hole-hole scattering contribution to the resistivity 
is given by 



V 



1 



\ 



(19) 



The latter matrix has a zero determinant and hence, in 
the absence of other scattering, allows for dissipationless 
flow of current of equal velocity in the two bands. 

We now focus on the dominant contribution, Eq. (|^). 
We begin with the calculation of the screened Coulomb 
interaction. The interaction and the response functions 
are described by 2 x 2 matrices, with indices correspond- 
ing to the two bands. Within the random phase approx- 
imation 



Vbarc 



where 



Vbaro — 



Vb Vb 
Vb Vb 



X 



bare 



Xu Xeh 

Xhe Xhh 



(20) 



(21) 



It follows from these equations that all screened in- 
teraction matrix elements are identical and given by 
Vsc{q,u}) = Vb{q)/e{q,uj), where Vb{q) = 2-Ee^leq is the 
bare interaction (e is the dielectric constant of the host 
material). The dielectric function is given by 



£(g,cj) = 
1 + [X«(9,w) 



Xihisi-,^) 



(22) 

XhM-, t^) + Xhh{q-, w)] V{q) . 



Plasmon excitations occur at wavevectors and frequen- 
cies which satisfy e[q,ujp{q) — ijp{q)] = 0, where ujp{q) is 
the plasmon frequency and jp^q) its width. To a good 
approximation, spin overlap may be ignored in the cal- 
culation of plasmons. Within this approximation, the 
intraband response functions are given by 



(18) X^^{q,^) 



(23) 



27r?i^ 



2Trh-' 



1 + i 



1 - 



(g,;«/c.)'-l 
1 - 



-1/2' 



-1/2 ~ 



(0 



for Lu < qv 



for LU > qVp^ 
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FIG. 23: Dashed arrow represents particle-hole excitation 
from the heavy to the light band at the minimal energy 
i^miniq) for a given momentum q (for q < k^p'' — k^P). This 
transition corresponds to excitation of a heavy hole at k'"P + q 
to the light band at k'"p^ in exactly the same direction in mo- 
mentum space. 



There are two plasmon branches in the system. The op- 
tical branch, with lj oc ,Jq, results from the solution of 

e{q,uj) = with uj > qvp^ in both bands. The acoustic 

(£) 

branch u q appears at frequencies, tu < qvp ■ 

Solving for the acoustic branch, we find to lowest order 
in ^ vnn /mh that £(q, bj) — for ijj(q) = Vpq — ijq, where 



Vp - i-f 



mh (h) 

Vp — I- 



2mi 



Ak 



(24) 



The screened interaction in the vicinity of the plasmon 
frequency is given by 



^. -• (25) 
^ — Vpq + «79 

Here, Vp is the interaction strength at the plasmon peak. 
Expanding the denominator to second order in a/ mf/m/j, 
we find 



ZTtTi Vp ' 



/ Wfe 

2mi 



(26) 



In the equations above we have used \Jmi/mh as a small 
parameter. At low densities this approximation breaks 
down but the plasmon dispersion remains linear and the 
interaction can still be approximated in the vicinity of the 
plasmon pole by eq. (|2^) with slightly different values 
of Vp^ 7, and Vp. The integral over frequency in (|^) is 
approximated in the vicinity of the plasmon frequency 
using 



frequency w may occur at wavevectors k satisfying Tiu 



|k-q 



. Possible k values are dictated by the S func- 



tion in the expression for Imxhe, 
ImXh?(q,c^) = ^ V(k. - q.)2|MM(k~ q,k)|^ 



k 



/o(^^„)-/o(4^01^( 



|k-q| 



(28) 



The possible phase space for interband particle-hole exci- 
tations is limited by a minimal frequency, uj^in{q). This 
frequency is illustrated in fig. ^ as the threshold for 
particle-hole excitation from the heavy to the light band. 
Such excitation is realized when the initial heavy hole 
momentum, k — q, and the final light hole momentum, 
k, are parallel and the final energy in the light band is 
equal to the Fermi energy (fig. |2^). The angle, 9, between 
k and q for such an excitation is vr. We are interested in 
excitations with frequencies close to ujmi^{q). Therefore, 
the angle 6 is close to tt, and A = 1 + cos is a small 
parameter. 

Using (|57]), the interband spinor overlap appearing in 
( p8| ) is given to leading order in A (and approximating 

k = k^p^) by 



|M;,,(k-q, k)|2 



rA, 



(29) 



where C, = ((d"^ |d''> + l/3(d"^ |d'^)) measures 

heavy/light hole mixing and its magnitude is approxi- 
mately 1 (exactly 1 in the absence of mixing). The inte- 
grals in (|2^) may now be calculated to leading order in 
A (using dO = d\/^/2X) to yield 



ImXM((7,cj) = 



3Ci 



and 



^v{q) = vi'^-v^;} 

kp +q 



nk^p\ 



All(30) 



(31) 



(32) 



Amax represents the maximal deviation of the particle- 
hole excitation from 6 = ir with k — kp. Note that 

liiiXhe{q,i^) oc [uj — Uminiq)]^^^ . For the calculation of 

the resistivity p\^^ we need Imx^|, eq. (|2^). Since the 



main contribution comes from a small range of k around 
—kp'^ci, we readily obtain 



\Vsciq,uj)\^ ^TT-/q\Vp\^S{u;-Vpq). (27) luYXhe{(l, uj) = [ky + q j cos^ iplmxhe{q,uj) , (33) 



We turn now to the calculation of the interband re- 
sponse function. An interband particle-hole excitation 
from the heavy to the light band at wavevector q and 



where ip is the angle between q and x. 

We now have all the ingredients needed for evaluat- 
ing Si, Eq. (|^). The interaction term, |Kicp7 is given 
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by Eq. (g^ ), so that all terms are calculated along the 
acoustic plasmon dispersion line, oj — Vpq. The interband 



response function Inix^^ is given by Eqs. (p 



_ 3|), hm- 

iting the integration range to q > qq (or u > ojq = VpQo). 
The intraband response of the light holes is given by 



Imx«(g,t^) 



mi 



(34) 



1 



j/qv 



Here, |M«|2 = |M«(k,k + q)|2 (||), calculated at a 
wavevector k contributing to Imx«(g,a;) (since huj <C 
Ep, the matrix elements are approximately the same). 
The integral (m) is calculated assuming T < fiujQ, hence. 



sinh"^(?iw/2r) w e 



/4. The integral is thus limited 



to the vicinity oi q = q^ and uj = loq. The temperature 
dependence of is readily evaluated. 



1 

T' 



dq {q - qo)3/2e-''''''(''-*)/^ 



90 



(35) 



!,k+q, £+1^ 



i",k"+q 



FIG. 24: Diagram for the Cooperon (see Eq. ^8|) The dashed 
line describes the short-range disorder. 



low temperatures. A full quantitative calculation takes 
all contributions into account. 

A contribution to the resistivity should arise from 
the different scattering processes at frequencies u < T. 
Preliminary calculations indicate this term should be ob- 
servable in the experiment. Its absence is hence surpris- 
ing. This discrepancy between theory and experiment 
might be related to one or more of the approximations 
involved in the calculation. These include the use of re- 
sponse functions calculated for q kp aX large wavevec- 
tors, q ~ 2kp, and the neglect of the anisotropy and 
nonparabolicity of the heavy band. 



We thus find that as a result of the minimal frequency 
for interband excitations, plj"* displays Arrhenius depen- 
dence upon T with a characteristic energy Tq — Tiloq. The 
magnitude and the T^/^-dependence of the prefactor are 
determined by details in the vicinity of the threshold for 
excitations. Collecting all terms and converting the scat- 
tering rate into resistivity one obtains 



(-1) 



1)^.V^ 



C 



(36) 



mi 



3/2 



5/2 



go 



T 
E, 



3/2 



The term C is given by 



1 



Ci|M, 



«| 



(37) 



l + i 



2m£ k 



(h) 



rrih k 



1 - 



V 



1 - 




We have considered here the plasmon enhanced con- 
tributio n to o ne of the scattering processes described in 
Section V A 2. With the help of Fig. pa we note that 



the only other plasmon enhanced scattering process in- 
volves carrier exchange between bands as described by 

(2) i k 

Pij , Eq. ( 1161) . The latter contribution comprises, how- 
ever, two terms that cancel each other within our approx- 
imation for the interband response function. Another 
possibility involving thermal activation of Im^ij beyond 
its zero temperature boundaries, in the Coulomb drag 
term, is discussed in Ref. ^ Thermal activation of plas- 
mons leads to an_Arrhenius dependence as well but with 
a T^/^ prefactor .E3 That contribution is hence smaller at 



B. Quantum interference corrections 

In this section we discuss quantum interference correc- 
tions to the conductivity. We assume a short-range inter- 
action and take the overlap between spin wavefunctions 
of the two bands, Eq. (57) into account. We find that in 
the limit of a small interband gap, the physics is similar to 
that of electrons in the presence of spin-orbit interaction. 
Quantum interference corrections i n this gfiS i g w ere pre- 
viously studied theoretically in refs.E30ca'EZlE3. Having 
calculated the Cooperon, the usual derivation (see, e.g, 
ref.E3) is followed to obtain the correction to the conduc- 
tivity. 

The following notation is used: Band indices of one 
particle are marked by Roman letter (i, j, . . .) that take 
the value £ or h. Greek letters {a,/3, . . .) denote band 
indices of the Cooperon two particles, taking the values 
£i, £h, h£ or hh. In this notation, the Cooperon is a 4 x 4 
matrix with Greek indices describing the bands of the 
incoming and outgoing particles. The Cooperon is given 
diagrammatically in fig. E4. At a total momentum q = 
and a total energy a; = it is given by (in this section 
we use units in which h = I) 

C„^(k,k") = V„^(k,k") (38) 
+ E / ^V„0(k,k')G+(k,e)G^:;(-k,e)C^^(k',k"), 



P 



The Green function is given by 



G±(k,e) 



.0) 



±i/2Tj- 



(39) 



where Tj is the mean free time in the j band, calculated 
from the disorder within the Born approximation. The 
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Green functions in Eq. ( p8| ) limit the contributing mo- 
menta to the vicinity of the Fermi surface. Taking the 
momenta on the Fermi surface, the integration over the 
magnitude k' is performed to give 



ith 



(40) 



/ (02G|+'(k,^)G^^(-k,6), (41) 



Finally, the disorder line is given by 

V„=,y, {9, 9') = {V^) Mr, {A9)Mj,, {A9) . (42) 

with (V^) being the disorder correlation function. Note 
the overlap between spin wavefunctions ( |57| ) is taken here 
on the Fermi surface and therefore depends on the angle 
difference A9 = 9' — 9 only. 

We discuss below two limits; A large band gap, £gT ^ 
1, as expected at high densities, and a small band gap, 
EgT <C 1, which becomes relevant as the density ap- 
proaches the metal-insulator transition. In both cases, 
for a short range disorder, the Cooperon equation can be 
solved exactly by calculating C — (l — Vn)~^V. 

We start with the high density limit, e^r ^ 1. In this 
case, the bare particle-particle propagator IIq, at q = 
0, w = is small for two particles in different bands. 
Neglecting this term, we find 



(2^ 



i^j. 

(43) 

The calculation reduces to 2 x 2 matrices with U and hh 
indices only. The mean free time is found in the Born 
approximation, using Eq. (|57|), to be 

^ = 2^^y' 0^G+ik',e){V')Mu'iA9)M,.i~A9) 

(^') /™ \^^<^)^ +(^v'M (44) 



The following observations simplify the calculation: In 
the equation for the Cooperon, (|o|), the disorder hue 
depends on A9 = 9' — 9, and 11^ is angle- independent. 
It is therefore possible to Fourier transform the equa- 
tion to angular momentum space. Using Eqs. ( |57| ) and 
(E3) , we note the only nonvanishing terms in Va0{A9) = 



T V 



"■^^ have < TO < 6. 



The Cooperon equation ( |40| ) takes the form, 

CW^VW+^V(">,C(™', (45) 



where Cap{9,9') = Em cS^e*"^^'"^^- The Cooperon 
has a large contribution whenever the matrix V^™^n^ 
has an eigenvalue close to 1. This happens for to = 3 
only, where 



/ 



V (d"'d" 



(46) 
J 



Since IIq — for a — £h or hi, we write V^^-* as a 2 x 

2 matrix with ii and hh indices only. Using (H), the 

(3) • • . — 

matrix V^^II^ has an eigenvalue 1, with eigenvector cx 

(-T-2,Tl). 

We show now that quantum interference corrections 
lead to weak antilocalization. Using the eigenvector of 

(3) 

V^^ n^, the Cooperon takes the form 



C{A9) 



T2/T1 -1 
-1 T1/T2 



J3A0 



[Dq^ - iLu)T 



(47) 



The diffusion constant D, and the mean free time r, re- 
flect the properties of both bands. In a weak localiza- 
tion diagram, each of the particle lines at one end of the 
Cooperon is connected with the other particle line at the 
other end of the Cooperon. Because the momenta of the 
two lines is opposite, = tt in the expression for the 
Cooperon and e^^^^ = — 1. This —1 term, which does not 
exist for the Cooperon in the single-band case, changes 
the overall sign of the diagram and leads to weak antilo- 
calization. 

In the limit of a large band gap we thus find weak an- 
tilocalization. We turn now to the opposite limit; a small 
gap, egT < 1. In that case we neglect differences in the 
density of states and scattering times in the two bands, 
as well as the bulk light hole contribution to the spin 
overlap functions, ^ The calculation is strictly a lead- 
ing order expansion in k. Here, the Cooperon equation 
( pq ) is not reduced to a 2 x 2 matrix since the propaga- 
tion of the two Cooperon particles in two different bands 
should be taken into account. The bare particle-particle 
propagator at q = a; = takes the form 



n = 



/I 




Vo 











\ 

' 

^ 

ie„T 

1/ 



(48) 



where indices are taken in the order i£, £h, hi, and hh. 
Although we assume CgT <gC 1, it is necessary to main- 
tain the egT terms in 11 as they lead to small deviations 
of the eigenvalues of VII from unity. Fourier decompo- 
sition in angular momentum space of the disorder line. 



(™)^imAe 





, has the following nonvanish- 
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ing components: 
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(49) 



(50) 



(51) 



We are interested in eigenvalues of IIV'^™-' which are close 
to unity and their corresponding eigenvectors. Calculat- 
ing C ~ {\ — Vn)^^V, we find the Cooperon 



(52) 




+ 



.1 -1\ 

{V^) [0 



-10 1 / 

/ 

1-10 

2 0-110 

V 



IU)]T 



„3iA0 




{Dq^ ~ iu + e'iT)T 



Evidently there are two spin-orbit scattering rates, t^J- = 
e^T and egT/2. When the Cooperon lines are matched in 
a weak localization diagram, A0 = tt. Consequently, a 

— 1 term appears in the second and third terms. The 

— 1 factor in the third term is cancelled against another 

— 1 factor appearing due to the contributing matrix el- 
ements, Cih^hi and C'hi,ih- The first, third and fourth 
terms hence generate weak localization while the second 
term contributes to weak antilocalization. The emerging 
picture is identical to the one obtained for spin-orbit scat- 
tering, where spin singlet states contribute to weak an- 
tilocalization while spin triplet states contribute to weak 
localization. 

The calculation of the quantum interference correction 
to the magnetoconductivity follows the standard route 



(see, e.g., Refs. leading to Eq. (|). This re- 

sult describes the crossover from weak antilocalization 



at egTTu 
in Fig. I 



> 1 to weak localization at e^TTip < 1, as seen 



Model and band structure 



In this section we discuss the band structure of a two- 
dimensional hole gas in a GaAs/AlGaAs heterostructuxe. 
This problem was considered theoretically in refs.oE3. 
Here, we present [the essential ingredients, relying on the 
derivation of ref.O. Our main purpose is to derive the 
overlap between spinor wavefunctions in the two bands, 
which is an important ingredient in the inter-band hole- 
hole scattering rate calculation. 

The lowest bulk valence bands of GaAs are fourfold de- 
generate at k = and split to light and heavy, doubly de- 
generate bands at finite k. When a confining potential in 
the z direction is introduced to create the 2DHG, the de- 
generacy between the bulk light and heavy holes is lifted 
at k = as well. In cases where the confining potential or 
the lattice lack inversion symmetry, spin-orbit coupling 
further lifts the twofold degeneracy of each band. The 
spin-orbit coupling depends on k thus at finite k, the va- 
lence band comprises four bands of different masses and 
Fermi velocities. Since the splitting between the bulk 
light and heavy holes is large, the relevant subbands for 
our experiment are the two originating from the bulk 
heavy band. We denote these bands as light and heavy 
but that nomenclature should not be confused with the 
light and heavy bulk bands. 

In the foUossfing calculation we use the soherical 
approximatioil£3. Consider the k-p HamiltoniarH within 
the spherical approximation, and in the presence of an 
asymmetric confining potential, 



H 



71 



57 

Y 



fc2 

2 



7(k- J)2 +l/(z). 



(53) 



Here, J are the angular momentum matrices for j = 3/2, 
and 71 = 6.85 and 7 = (272 -I- 373)/5 = 2.58 are the Lut- 
tinger parameters for GaAs in atomic units. The confin- 
ing potential V{z) is usually caleulated self-consistently 
within a Hartree approximationoE3 . Here it is sufficient 
to assume that V{z) is the result of such a calculation. 

A note regarding our notation: The bulk light and 
heavy holes are denoted by capital letters {L and H). 
The light and heavy bands originating from the bulk 
heavy band are denote by £ and h. It is the latter two 
which are relevant for our discussion. 

Within the subspace of light and heavy bulk states, 
and using a J = | represe ntati on the two dimensional 
Hamiltonian takes the formEJEj, 
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(54) 



Here, fc± = kx^iky. The parameters s = ('0t*'''|V'ff'''): a-nd 



i/C = (^'1°'' Ifczl^/'^''), where iPh't'^l'' are the z-direction 
states corresponding to the bulk heavy and Hght holes, 
respectively reflect the wave function in the third dimen- 
sion. Note K ^ for asymmetric confining potential. 
To third order in k the Hamiltonian's eigenvalues are 



/,(o)\ 



/,(o) JO) 



71+7 



.(0) 



k'± 



M JO)' 



(55) 

with the upper (lower) sign corresponding to i = £ {i = 
h). Within this model the two bands are degenerate in 
the case of a symmetric confining potential {K = 0). 

The eigenvectors of iJk corresponding to the two rele- 
vant bands are U{0)'E!^^'^ and U{6)'E^'^\ where 



u{e) 



1 







\ —ie 



3ie 



1 \ 



3ie 



le 



(56) 



the angle 9 marks the direction of k in the x- 
y plane (this transformation is a single valued ver- 
sion of the unitary transformation used in ref. 34), 



-e^^«d'^(fc),0,0 



and S('')(k) = 



sW(k) = (d^)(fc), 

(o,0,e'^'''d^\k),^i!;Hk)y Here are chosen to be 

real, and the dependence upon the direction of k is given 
explicitly in the expressions for 2'^*^(k). The ^'s, depend- 
ing only on the magnitude of k, may be explicitly calcu- 
lated, but are not needed here. For our calculation we 
need the overlap of the spin part of two different eigen- 
states of the system, as a function of the two momenta 
directions. This overlap is given by 



My(k,k') = ^^<''\k)U\e)U{9')^<^^\k') 



(57) 
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FIG. 25; f{Si,Sh,Q,H) vs. magnetic field at two different 
temperatures for ptotai = 4 x 10^^ cm~^ . 



where A9 9' — 9, and the indices i and j take the 
value £ or h. This result is used in the calculation of the 



transport properties in Section (VA2) 
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APPENDIX A 

The function f{Si,Sh, Q, H) in Eq. (|l|) is given by 
1 1 



[1 + [H/Wr]^ (Si + Sh + 2Q)4 
[Co + Ci-H^ + C2-H^ + C3- i?'] , 

where 

Co = -iSi + Sh + 2Qf -{SiSh-Q^, 
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Ci = Qo + Q-Qi + Q^ ■Q2 + Q^ -Qs + Q^- Qa, Qz = 8{RlSi + RfSn), 
C2 = Po + Q-Pi + Q^- P2, Qi^ GRiRh + Rl, 

C3 — Rf Rf^{Ri + Rh)^ , Po — 2RfRj^{Si ~ Sh)^ ~ [R^Si — RfSh)^ 

+ 2RiRh{RlSf + RfSl), 
Pi = ARiRh[RhSi{Ri + 2Rh) + RiSh{Rh + 2i?0], 
Qo = (RhSf + RiSlf ~ 2{Si + Sh){RlSi + RfSh)SiSh, P^ = 2RiRh{RI + iif). 
Qi = A[{Si~Sh){RlSt-R^iSl) + 2RiRhSiSh{Si + Sh)\, 

Q2 = 2 [SiSh{Ri + Rhf + RiRh{Si - ShY -p.g n depicts f{Su Sh, Q, H) vs. magnetic field at T : 

+ A{RhSi + RiShf] , 0.4 and 1.5K for ptotai = 4 x 10"cm-2. 
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